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Abstract 

We study B K(*k+e- decays {£ = e, n) based on a data sample of 657 million BB pairs 
collected with the Belle detector at the KEKB e"'"e~ collider. We report the differential branching 
fraction, isospin asymmetry, K* polarization, and the forward-backward asymmetry {Afb) as 
functions of = Mff^cP'. The fitted Aps spectrum exceeds the Standard Model expectation by 
2.7 standard deviations. The measured branching fractions are B{B K*i'^£~) = {10.7~^\'q ± 
0.9) X 10~'^ and B{B Ki+i') = (4.8lo:4 ± 0.3) x 10"^ where the first errors are statistical and 
the second are systematic, with the muon to electron ratios Rk* = 0.83ib0.17 it 0.08 and Rk = 
1.03±0.19 ± 0.06, respectively 

PACS numbers: 13.25 Hw, 13.20 He 
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The b si^i~ transition is a flavor-changing neutral current (FCNC) process, which, 
in the Standard Model (SM), proceeds at lowest order via either a Z/7 penguin diagram 
or a W^W~ box diagram. Since their amplitudes may interfere with the contributions 
from non-SM particles [l|, the transition can probe the presence of yet unobserved particles 
and processes. More specifically, the lepton forward-backward asymmetry {Aps) and the 
differential branching fraction as functions of dilepton invariant mass {Mjp) in the decays 
B — > K*i~^i~ differ from the SM expectations in various extended models [2(]. The former is 
largely insensitive to the theoretical uncertainties of the form factors describing the decay, 
and can hence provide a stringent experimental test of the SM. The latter has been so far 
determined only with a modest precision 0, 0] • It can be used to extract the information 
on the coefficients associated with the theoretical models as well. 

In this paper, we report measurements of the differential branching fractions and the 
isospin asymmetries as functions of = Mf^c^ for B — > K*i~^i~ and B Ki~^i~ decays. 
The K* polarization and ApB for B K*i~^i~ decays as functions of are presented as 
well. A data sample of 657 million BB pairs, corresponding to 605 fb~^, collected with the 
Belle detector at the KEKB asymmetric-energy e~^e~ collider is examined. The Belle 
detector is described in detail elsewhere j6|. 

We reconstruct B — >■ K^*H'^£~ signal events in 10 final states: K^7i~, i^^vr"*", K~^7i^, 

combined with either electron or muon pairs. All charged tracks 
other than the — > vr+vr" daughters are required to be associated with the interaction 
point (IP). A track is identified as a (tt"*") by combining information from the aerogel 
Cherenkov and time-of-fiight subsystems with dE/dx measurements in the central drift 
chamber [8|. The kaon (pion) identification is more than 85% (89%) efficient while removing 
more than 92% (91%) of pions (kaons). For muon (electron) candidates, we use the lepton 
identification likelihood described in Ref. jl], which we denote by TZx [x denotes or e). We 
select /i^ candidates with > 0.7 GeV/c and momentum-dependent 7^^ requirements that 
retain (93.4%±2.0%) of muons while removing (98.8%±0.2%) of pions. Electron candidates 
are required to have Tie > 0.9, 71^ < 0.8, and pe > 0.4 GeV/c. These requirements retain 
(92.3%±1.7%) of electrons while removing (99.7%±0.1%) of pions. Bremsstrahlung photons 
emitted by electrons are recovered by adding neutral clusters found within a 50 mrad cone 
along the electron direction. The cluster energies are required to be between 20 and 500 MeV. 

Pairs of oppositely-charged tracks are used to reconstruct Kg 'iT'^7r~ candidates. 
The invariant mass is required to lie within the range 483-513 MeV/c^ (±5 times the Kg 
reconstructed- mass resolution). Other selection criteria are based on the distance and the 
direction of the Kg vertex and the distance of daughter tracks to the IP. For 7r° — s> 77 
candidates, a minimum photon energy of 50 MeV is required and the invariant mass must 
be in the range 115 MeV/c^ < M^j < 152 MeV/c^ (±3 times the 7r° reconstructed-mass 
resolution). Requirements on the photon energy asymmetry \E}^ - E^\/{E}^ + E^) < 0.9 and 
vr^ momentum p^o > 200 MeV/c suppress the combinatorial background. 

S-meson candidates are reconstructed by combining a K^*^ candidate and a pair of 
o ppositely cha rged leptons, and selected using the beam-energy constrained mass Mbc = 
V -^beam ~ Pb ^^"^ ^he energy difference AE = Eb — -Ebeam, where Eb and pb are the re- 
constructed energy and momentum of the B candidate in the T(4S') rest frame and -Ebeam 
is the beam energy in this frame. Bremsstrahlung photons are included in the calculation 
of the momenta of electrons and hence are included in the calculations of Mbc, -^'beam and 
g^. We require i?-meson candidates to be within the region Mbc > 5.20 GeV/c^ and —35 
(—55) MeV < AE < 35 MeV for the muon (electron) modes. The signal region is defined as 
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5.27 GeV/c^ < Mbc < 5.29 GeV/c^. For the K* modes, the Mk-k candidate (signal) region 
is defined as Mk^ < 1.2 GeV/c^ (|Mx^ - m^*! < 80 MeV/c^). 

The main backgrounds are continuum e~^e^ i<l = d, c, s) and semileptonic B 

decay events. We use the same set of variables and the likelihood ratio, TZ, described in 
Ref. [9I, for continuum e'^e~ qq suppression. For the suppression of semileptonic B 



decays, we combine a Fisher discriminant including 16 modified Fox- Wolfram moments [10 



the missing mass Mmiss = \/ -E'^iss ~ Pmiss' cosOb, and the lepton separation near the IP in 
the z direction to form the likelihood ratio T^si = 'Cs/(£s + -^^si), where -Emiss(Pmiss) is the 
missing energy (momentum), 6b is the polar angle between the reconstructed B candidate 
and the beam direction in the T{4S) rest frame, and Cg {Cs\) is the likelihood for signal 
(semileptonic B) decays. Combinatorial background suppression is improved by including 
q^ and 5-fiavor tagging information jll|. Selection criteria for TZ and T^si are determined 



by maximizing the value of S/ y/S + B, where S and B denote the expected yields of signal 
and background events in the signal region, respectively, in different and tagging regions. 

The dominant backgrounds that peak in the signal region are from B ^ J/ ipX and ip'X 
and rejected in the following q^ regions (in units of GeV^/c^): 8.68 < q'^{fi~^fi~) < 10.09, 
12.86 < g2(/x+/i") < 14.18, 8.11 < q^{e+e-) < 10.03, and 12.15 < q^{e+e-) < 14.11. To 
remove the background from B J/iIj{iIj')K* events with one of the muons misidentified as a 
pion candidate, we reject events with —0.10 GeV/c^ < M^^ — m < 0.08 GeV/c^, where 
the pion is assigned the muon mass. Background from B —>■ DX is rejected by additional 
veto windows |Mx^ — m£)| < 0.02 GeV/c^ and \MKT,^ — mB,\ < 0.02 GeV/c^, where the muon 
is assigned the pion mass. The invariant mass of an electron pair must exceed 0.14 GeV/c^ 
in order to remove background from photon conversions and 7r° 7e^e~ decays. 

If multiple B candidates survive these selections in an event, we select the one with the 
smallest jAE'l. The fractions of multiple B events are about 7%, 12%, and 20% for the 
K^n~ , K^iT~^, and K~^tt^ modes, respectively, and less than 1% for the and Kg modes, 
according to a study using Monte Carlo (MC) simulation of 5 — > K^*H£ decays. 

To determine the signal yields, we perform an extended unbinned maximum likelihood 
fit to Mbc and Mk^ for B K*i+i- decays, and to Mbc for B Ki+i- decays. The 
likelihood function includes contributions from signal, combinatorial, B — > J/ilj{ip')X , and 
B K^*\7r backgrounds. The signal PDFs consist of a Gaussian (Crystal Ball func- 



tion [12[) in Mbc for the muon (electron) modes and a relativistic Breit-Wigner shape in 



Mktt for the K* resonance. The means and widths are determined from MC simulation and 
calibrated using B — *• J/tpK^*^ decays. The PDFs of signal decays in which either the kaon 
or the pion candidate is wrongly associated to the K* decay (self-cross-feed) are modelled 
by a two-dimensional smoothed histogram function with g^-dependent fractions obtained 
from MC simulation of the signal decays. The combinatorial PDFs are represented by the 
product of an empirical ARGUS function [l3| in Mbc with the sum in M^n of a threshold 
function (whose threshold is fixed at rriK + "^7r) and a relativistic Breit-Wigner shape at the 
K* resonance. The PDFs and yields for B J/i(j{i(j')X decays are determined from a large 
MC sample, while the B — » K^*^7nr PDFs and normalizations are determined from data, 
taking into account the probabilities of pions to be misidentified as muons. Yields for signal 
and combinatorial background and the combinatorial PDF parameters are allowed to fioat 
in the fit, while the yields and parameters for other components are fixed. 

We divide q^ into 6 bins and extract the signal and combinatorial background yields in 
each bin. The K* longitudinal polarization fractions (Fl) and ApB are extracted from fits 
to cos6'x* and cos 9b£, respectively, in the signal region, where 9k* is the angle between the 
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kaon direction and the direction opposite to the B meson in the K* rest frame, and Ose is 
the angle between the and the opposite of the B (B) direction in the dilepton rest 

frame. The signal PDFs for the fit to cos 6k* and cosObi are described by a product of the 
i^T* /dilepton polarization function and the efficiency, 

[fFiCos^^x. + |(l-FL)(l-cos2^i^0 ] xe^cosOK*) 

and 

[ - cos2 Ob^) + 1(1 - + cos2 OBd 

+Afb cos 9b£ ] X e{cos 9bi) , 

respectively. The first two terms in the dilepton polarization function correspond to the pro- 
duction of K*^s with longitudinal and transverse polarization, while the third term generates 
the forward-backward asymmetry. Figures in Ref. [13] illustrate the fits for B yields, Fl, 
and ApB in each bin. In the fit to cos 9k* {cos9Be), Fl {Afb) is the only free parameter, 
while the other PDFs and normalizations are fixed. For the B Ki~^i~ modes, we set 
Fl = 1 and the B Kgi~^i~ sample is not used. 



TABLE I: Fit results in each of six q'^ bins and an additional bin from 1 to 6 GeV^/c^ for which 
recent theory predictions are available 15|]. The first uncertainties are statistical and the second 
are systematic. 
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Table [T] lists the measurements of B yields; the partial branching fractions, obtained by 
correcting the B yields for g^-dependent efficiencies; Fl; and Afb individual bins. In 
the calculation of the partial branching fractions, we adopt the SM lepton fiavor ratios of 
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the muon to electron modes 1^ and express the branching fractions in terms of the muon 
channel. The ratio for the full interval is = 0.75 = 1) for the B K*e+i- 

{B — > Ki~^i~) mode, where the deviation from unity is due to the photon pole. The results, 
as well as the SM curves, are shown in Fig. [H To illustrate how non-SM physics might 
manifest itself, we superimpose curves on the and ApB plots corresponding to the case 
of Cy with reversed sign (C7 = —Cj^^). The measured values do not reject this possibility. 

The total branching fractions, extrapolated from the partial branching fractions 
measured to be 



are 
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B{B - 
B{B 



ICtt 
^ Ktt 



(10.7i};J±0.9) X 10"^ 
(4.8lJ]i ± 0.3) X 10"^ ; 



while the fitted CP asymmetries, defined in terms of the B (B) yield A'^^ (A'^^) as Acp = 
im-N-,)/iNt + Nj,), are Q 



AcpiK* 
Acp{K+tl 



-0.10 ±0.10 ±0.01 
0.04 ±0.10 ±0.02 . 



The lepton flavor ratio is sensitive to Higgs emission and is predicted to be larger than 
the SM value in the Higgs doublet model at large tan/3 [l8|. The measured ratios are 

Rk* 
Rk 



0.83 ±0.17 ±0.08 , 
1.03 ±0.19 ±0.06 . 



The isospin asymmetry, shown in Table [T] and Fig. [H is defined as 



Ai 



{tb+Ztbo) X B{K(*)0£^ 



B{K(*)-- 



where tb+/tbo = 1.071 is the lifetime ratio of B^ to B^ [l9[. A large isospin asymmetry for 
below the mass of the J/^ resonance was reported recently [20|. We also measure the 
combined Aj for < 8.68 GeV^/c^ and find 



Aj{B K*tt 
Ai{B Ktt 
Ai{B K^*k+t 



-0.29l°:J^ ± 0.09 



1+0.17 
'--0.14 

1+0.12 



-0.31^^;i^±0.08 



cr 



a 



0.30I^;iJ ± 0.08 a 



1.37 
1.75 
2.22 



where a denotes the significance from a null asymmetry and is defined as a = 
y^— 21n (/^o/'^max), where £0 is the likelihood with Aj constrained to be zero and £max is the 
maximum likelihood. Systematic uncertainties are considered in the significance calculation. 
No significant isospin asymmetry is found at low q^. 

Systematic uncertainties in the branching fraction measurements arise predominantly 
from tracking efficiencies (2.0%-4.4%), MC decay models (0.9%-4.6%), electron (3.0%) and 
muon (2.6%) identification, (4.9%) and vr" (4.0%) reconstruction, and 71 and TZ^i selection 
(1.2%-3.6%). The MC simulated samples of the signal are generated based on a model 
derived from Ref . 1^ . The modeling uncertainties are evaluated by comparing MC samples 
based on different models |2l|, while lepton identification is studied using a. J/ip —>■ i~^£~ data 
sample. For TZ and T^si selections, we estimate the uncertainties from large control samples 
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FIG. 1: Differential branching fractions for the (a) K*i^£~ and (b) Ki'^i^ modes as a function 
of g^. The two shaded regions are veto windows to reject J/'ijj{'ip')X events. The sohd curves show 
the SM theoretical predictions with the minimum and maximum allowed form factors [la], (c) and 
(d) show the fit results for and Aps in K*i'^£~ as a function of q"^, together with the solid 
(dotted) curve representing the SM (C7 = —C^^) prediction 16[]. (e) is the Aj asymmetry as a 
function of for the K*£^£~ (filled circles) and K£~^£~ (open circles) modes. 

with the same final states, B J/ipK^*"* with J/ip ^ Other uncertainties such as 

kaon and pion identification efficiencies, fitting PDFs, background contamination from J/ ijj 
decays and charmless B decays, and the number of BB pairs are found to be small. The 
total systematic uncertainties in the branching fractions for different decay channels are 
6.8%-12.2% and 5.2%-7.4% for the K*i+i- and Ki+i- modes, respectively. 

The main uncertainties for angular fits are propagated from the errors on the fixed nor- 
malizations and Fl, determined from M]^c^Mxn and cosOk* fits, respectively. Fitting bias 
and fitting PDFs are checked using large B J/ipK^*^ and MC samples. The total un- 
certainties for the Fl and Aps fits depend on the bin and range from 0.02-0.06 and 
0.03-0.13, respectively. The systematic errors on Acp are assigned using the CP asymme- 
try measured in sideband data without TZ and T^si selections and are found to be 0.01-0.02. 
The systematic error on -R/^(.) {Aj) is determined by combining the uncertainties from lep- 
ton (K/tt) identification, TZ and Tls\ selections, fitting PDFs and background contamination. 
The uncertainty in Aj from the assumption of equal production of B^B^ and B^B~ pairs 
is also considered. The correlated systematic errors among bins are negligible for all the 
measurements. 

In summary, we report the differential branching fraction, isospin asymmetry, K* longi- 
tudinal polarization and forward-backward asymmetry as functions of g^, as well as total 
branching fractions, lepton fiavor ratios, and CP asymmetries for B K^*'>£~^£~. These 
results supersede our previous measurements jl] and are consistent with the latest BaBar 



results 1^, |20| with better precision. The differential branching fraction, lepton flavor ratios, 
and K* polarization are consistent with the SM predictions. No signiflcant CP asymmetry 
is found in the study. The isospin asymmetry does not deviate signiflcantly from the null 
value. The Apsiq"^) spectrum for B K*£^£~ decays tends to be shifted toward the pos- 
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itive side from the SM expectation. The cumulative difference between the SM prediction 
and the measured points is found to be 2.7 standard deviations. A much larger data set, 



which will be available at the proposed super B factory [22|] and LHCb |23|, is needed to 
make more precise tests of the SM. 
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APPENDIX 



Below is EPAPS supplementary material for fit figures, as well as the total branching frac- 
tions and CP asymmetries for individual K^*H~^i~ modes, which is mentioned in Ref. 14 



TABLE II: Total branching fractions and CP asymmetries for subsamples of B 
B Ki+i~ decays. 



K*£+i~ and 



Mode 


B (10-7) Acp 




ll.l+l^ibl.O -0.12+°;^^±0.02 
10.6li;^±0.7 0.00l[j:J^±0.03 

ii.oli:^±o.8 -o.o3lg:J^±o.o2 


K*+e+e~ 
K*e+e- 


i7.3j:4:2±2.o -o.utoil±om 

ll.Sll^zbO.g -0.21+°;J^±0.02 
13.9l|g±1.2 -0.18lg:Jf±0.01 




12.4t|?±1.3 -0.13toil±Om 
9.7l|;5±0.7 -0.08l[]:J^±0.02 

io.7ti:o±o.9 -o.ioj:[!:J[]±o.oi 




5.3+[5:^±0.3 -0.05t[!:Ji^±0.03 



4.4+|:?±0.3 
5.0l[]:^±0.3 



K+e+e- 



5.7lJ]|±0.3 
2.0l|:^±0.1 
4.8+55:^±0.3 



0.14to:M±0.03 



K£+ 



5.3lJ]|±0.3 
3.4l°|±0.2 



o.o4j:o:}o±o.o2 
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FIG. 2: Fits to Mktt, ^h^, cos Bk*, and cos 6*5^ for the K*i+£~ decays in 6 q"^ (GeVVc^) bins: 
(a)~(d) 0.00-2.00, (e)~(h) 2.00-4.30, (i)~(l) 4.30-8.68, (m)~(p) 10.09-12.86, (q)~(t) 14.18-16.00, 
and (u)~(x) > 16.00. The sohd, long-dashed, glprt-dashed, and dotted curves represent the com- 
bined fit result, fitted signal, combinatorial background, and J/il:[Tp')X background, respectively. 




M^^ (GeV/c^) M^^ (GeV/c^) 



FIG. 3: Fits to Afbc and cosO bi for the Kl+i- decays in 6 (GeVVc^) bins: (a)~(b) 0.00- 
2.00, (c)~(d) 2.00-4.30, (e)~(f) 4.30-8.68, (g)~(h) 10.09-12.86, (i)~(j) 14.18-16.00, and (k)~(l) > 
16.00. The soHd, long-dashed, and short-dashed represent the combined fit result, fitted signal, and 
combinatorial background, respectively. Th.eJ /ip{ip')X background is so small and not represented. 
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